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ABSTRACT 


This  report  presents  an  algorithm  and  HP-67/97  programs 
for  position  determination  with  Loran-C  chains.  Operational 
data  cards  are  prepared  in  advance  for  Loran-C  triplets. 
Position  determination  is  performed  using  a  single  program  card 
and  an  appropriate  operational  data  card. 


A.  Introduction 


The  Loran  system  is  a  radio  aid  to  navigation  which 
utilizes  the  principle  of  hyperbolic  fixing.  The  locus  of 
points  for  which  the  difference  in  arrival  time  of  synchronized 
signals  from  a  pair  of  transmitters  is  constant  determines  a 
hyperbolic  line  of  positions  (LOP) .  The  intersection  of  two 
hyperbolic  lines  of  position  from  two  pairs  of  transmitters 
determines  position  or  a  hyperbolic  fix.  That  two  pairs  of 
stations  are  required  for  a  fix  does  not  necessarily  mean  that 
there  are  four  separate  stations,  for  one  station  of  one  pair 
may  be  colocated  with  one  station  of  the  other  pair  forming 
a  Loran  triplet  (Figure  1) .  Triplets  may  be  joined  "end-to-end" 
by  station  colocation  to  form  a  Loran  chain  (Figure  2) .  Loran 
chains  are  common  on  both  the  East  and  West  Coasts  of  the 
North  American  continent. 

The  early  "Standard  Loran"  or  Loran-A"  operating  at  a 
frequency  just  below  2MHz  is  still  in  use  in  the  Pacific  area. 
The  present  day  "Loran-C"  operates  at  100-kHz  and  is  in  use 
in  both  the  Atlantic  and  Pacific  Areas.  The  computational 
algorithm  and  programs  described  herein  can  be  used  for  position 
determination  with  Loran-C  triplets.  Further  information  on 
the  history,  development  and  operation  of  the  Loran  systems 
may  be  found  in  References  1  and  2. 


(a)  Colocated  Master  Stations 


(b)  Colocated  Slave  Stations 


(c)  Colocated  Master  and  Slave 


Figure  1.  Loran  Triplets 


B.  Program  Description 

One  program  card  and  one  operational  data  card 
(described  below)  are  all  that  is  required  for  on-location 
position  determination  from  Loran  triplet  time-difference 
measurements.  Two  program  cards  are  required  to  prepare 
operational  data  cards;  these  operational  data  cards  should 
be  prepared  and  validated  prior  to  on-location  navigational  use. 
Thus  although  three  program  cards  are  described  only  one  program 
card  is  required  for  navigation;  two  program  cards  are  used  to 
prepare  operational  data  cards  during  or  prior  to  mission  planning. 
The  function  of  each  program  card  and  its  intended  use  follows. 


Program  Card  1 .  This  program  card  is  used  to  prepare 
master >  data  cards.  A  master  data  card  requires  the  following 
information  for  a  master (M)  station/slave (S)  station  pair: 

1.  A  M/S  pair  identification  number. 

2.  The  quantity  At  which  is  the  sum  of  the  coding  delay  plus 
the  one  way  base  line  time  in  microseconds. 

3.  The  latitude  and  longitude  of  the  master  station. 

4.  The  latitude  and  longitude  of  the  slave  station. 

Some  preprocessing  of  these  data  is  performed  before  the  master 
data  card  is  generated.  The  data  generated  require  only  one  side 
of  an  HP-67/97  magnetic  card  for  each  M/S  pair,  thus  a  second 
M/S  pair  may  be  placed  on  side  2  of  the  card  (thus  conserving 
cards)  if  desired.  It  is  envisaged  that  a  master  data  card 
will  be  prepared  in  advance  for  each  M/S  pair  that  might  be 
received  within  an  area  of  operation. 
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Program  Card  2 .  This  program  card  is  used  to  prepare 
an  operational  data  card  for  every  Loran  triplet  within  an 
operational  area.  Each  operational  data  card  contains  data 
merged  from  the  master  data  cards  which  contain  M/S  pair 
information  for  each  pair  of  the  triplet.  These  merged  data 
are  validity  checked,  colocation  of  master  or  slave  determined 
and  encoded. 

The  only  inputs  required  for  this  program  are  the  two 
master  data  cards  that  comprise  the  Loran  triplet.  It  is 
possible  to  prepare  and  store  operational  data  cards  rather 
than  master  data  cards.  This  may  be  desirable  if  there  is  no 
scarcity  of  cards  and  storage  space,  however  the  number  of 
possible  Loran  triplets  is  considerably  larger  than  the 
number  of  M/S  pairs. 

Program  Card  3.  This  program  card  is  used  in  conjunc¬ 
tion  with  an  operational  data  card  for  position  determination. 
Required  input  is  the  indicated  time  difference  T  for  each 
M/S  pair  of  the  triplet.  Output  is  the  computed  latitude  and 
longitude  of  the  fix.  Note:  Every  Loran  fix  has  two  possible 
solutions.  The  unwanted  solution  can  almost  always  be  rejected 
by  inspection,  however,  if  the  stations  of  the  Loran  triplet 
are  nearly  aligned  then  either  solution  may  be  valid  even 
though  only  one  solution  should  be  consistent  with  the  flight 
plan . 
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HP-67/97  Calculator  Programs 
1.  User  Instructions 


CARD  1 


Input 

1 

Step 

Instructions 

Data/Units 

Keys 

1. 

Read  in  program  card  (both 
sides) 

2. 

Input  a  unique  ID  number  for 
the  Loran  pair* 

ID 

f  a 

3. 

Input  the  coding  delay  At 

At 

f  c 

4a. 

Input  the  master  station  latitude 

4>m 

+ 

b. 

and  longitude  (CHS  for  West)** 

xm 

A 

5a. 

Input  the  slave  station  latitude 

4>s 

t 

b. 

and  longitude  (CHS  for  West) 

xs 

C 

6a . 

Run 

None 

E 

b. 

Pass  a  blank  data  card  through 
the  card  reader. 

Output 

Data/Units 


Note :  Loran  pairs  are  coded  on  the  navigation  maps  using 

designators  such  as  9930X,  9930Y,  99303  and  9930W.  It  is 
suggested  that  the  ID's  for  these  pairs  be  coded  as 
9930.1,  9930.2,  9930.3  and  9930.4,  respectively.  However, 
any  consistent  scheme  is  acceptable. 


The  format  for  position  data  input  is  of  the  form: 
+DDD.MMSSFF,  where 

ODD  denotes  degrees 
MMM  denotes  minutes 
SS  denotes  seconds 
FF  dnotes  hundredths  of  a  second. 

The  minus  sign  (-)  denotes  Southern  latitudes  or  Western 
longitudes . 


CARD  2 


I 

i 


Step 

Instructions 

Input 

Data/Units 

Keys 

Output 

Data/Units 

1. 

Read  in  program  card  2 

2a  • 

Start 

— 

E 

9.00 

b. 

9.00  will  flash  in  the  display. 
Insert  a  master  data  card 
containing  the  first  pair  of 
a  Loran  triplet. 

c. 

9.00  will  flash  in  the  display 
once  more.  Insert  a  master 
data  card  containing  the  second 
pair  of  the  Loran  triplet. 

9.00 

d. 

If  the  data  form  a  proper 
triplet,  "crd"  will  appear  in 
the  display. 

crd 

e. 

Pass  both  sides  of  a  blank  card 
thru  the  card  reader  to  produce 
the  operational  data  card  for 
the  Loran  triplet. 

Should  "error"  appear  in  the  display,  then  the  two  master  data 
cards  do  not  compare  to  form  a  Loran  triplet.  Both  the  latitude 
and  longitude  of  the  colocated  stations  must  be  identical  on 
both  master  data  cards  in  order  to  successfully  produce  an 
operational  data  card. 


Label  the  A  key  position  with  the  identification  number  of  the 
first  Loran  pair  (the  pair  inserted  in  Step  2b)  and  label  the  B 
key  position  with  the  identification  number  of  the  second  Loran 
pair  (from  Step  2c). 


CARD  3 


Step 

Instructions 

1. 

Read  in  both  sides  of  the  program 
card  3. 

2. 

Read  both  sides  of  the  operational 
data  card  for  the  Loran  triplet 
that  you  are  receiving. 

3a. 

j  Set  to  compute  Solution  A. 

b. 

Set  to  compute  Solution  B. 

4. 

Input  the  observed  time  delay 
from  the  first  Loran  pair. 

5. 

Input  the  observed  time  delay 
from  the  second  Loran  pair. 

6. 

|  Compute  fix  Latitude 

Longitude 

l 

7. 

J  Repeat  from  Step  2  with  a  new 

1  operational  data  card  or  from 

Steps  3  or  4  as  required. 

Input 

Data/Units 


Output 

Keys  Data/Units 


Latitude 

Longitude 


M 


m 

Input 

Output 

Step 

Instructions 

Data/Units 

Keys 

Data/Units 

In  this  series  of  examples  we 
will  prepare  and  use  data 
cards  for  the  Loran-C  pairs 
9930X  and  9930Y. 

1. 

Read  in  program  card  1 
(both  sides) 

2. 

Input  the  ID  for  9930X 

9930.1 

f 

a 

9930.10 

3. 

Input  the  coding  delay  At 
for  9930X 

36389.66 

f 

c 

36389.66 

4a. 

Input  the  master  station 
latitude 

34.034604 

+ 

b. 

and  longitude  (CHS  for  West) 

-77.544676 

A 

— 

5a. 

Input  the  slave  station 
latitude 

46.463218 

1 

b. 

and  longitude  (CHS  for  West) 

-53.102816 

C 

— 

6a. 

Compute 

None 

E 

crd 

b. 

Pass  a  blank  data  card  through 
the  card  reader.  Label  the 
card  9930X  MASTER 

— 

7. 

Input  the  ID  for  9930Y 

9930.2 

f 

a 

9930.2 

8. 

Input  the  coding  delay  At 
for  9930Y 

52541.31 

f 

c 

52541.31 

9a. 

Input  the  master  station 
latitude 

34.034604 

+ 

b. 

and  longitude  (CHS  for  West) 

-77.544676 

A 

— 

10a. 

Input  the  slave  station 
latitude 

41.151193 

t 

b. 

and  longitude  (CHS  for  West) 

-69.583909 

C 

— 

11a. 

Compute 

None 

E 

crd 

b. 

Pass  a  blank  data  card  (or 
the  second  side  of  the  card 
used  in  Step  6b)  through  the 
card  reader.  Label  the  side 
9930Y  MASTER 

12. 


These  t-wn  cards  will  be  used 
in  the  next  example. 


CARD  2 


Step 

Instructions 

— 

Keys 

Output 

Data/Units 

1. 

Reed  in  program  card  2. 

2a. 

Start 

— 

E 

9.00 

b. 

While  9.00  is  flashing  in  the 
display,  insert  the  MASTER 
data  card  for  station  9930X 
into  the  card  reader. 

c. 

When  9.00  starts  flashing  in 
the  display  again,  insert 
the  MASTER  data  card  for 
station  9930Y  into  the  card 
reader. 

9.00 

d. 

"Crd1  will  appear  in  the 
display. 

— 

crd 

e . 

Pass  both  sides  of  a  blank 
card  through  the  card 
reader. 

Label  this  card  9930X/9930Y 
OPERATIONAL  DATA  CARD.  Then 
label  the  A  key  position  9990X 
and  the  B  key  position  9930Y. 
This  card  will  be  used  in  the 
next  example. 

0.00 
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CARD  3 


Step 

Instructions 

Input 

Data/Units 

Keys 

■ 

You  are  receiving  9930X  and 

9930Y  and  wish  to  obtain  a 
fix. 

B 

Read  in  program  card  3 
(both  sides) 

B 

Read  in  the  operational  data 
card  for  the  triplet 

9930X/9930Y  (both  sides) 

Bl  B 

Select  Solution  A. 

f  a 

B 

The  indicated  time  delay  is 

49400  Ms  from  9930Y. 

Input  the  indicated  time  delay. 

49400 

B 

5. 

The  indicated  time  delay  is 

28800  ms  from  9930X. 

Input  the  indicated  time  delay. 

28800 

A 

6. 

Solution  A:  42 ° 44 ' S7"N  Latitude 
41o07'32"W  Longitude 

E 

R/S 

Solution  B:  27°00' 07"s  Latitude 
102°27'12"E  Longitude 
Since  you  are  navigating  over 
the  North  Atlantic,  Solution  A 
is  the  desired  fix. 

■ 

Repeat  from  Step  2  with  a  new 
operational  data  card  or  from 
Steps  3  or  4  as  required. 

0.00 

0.00 

42.4457 

-41.0732 
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3.  Program  Storage  Allocations  and  Program  Listings 


Card  1. 

Registers: 

RO: 

ID 

SO: 

L 

RA: 

Rl: 

At 

SI: 

T 

RB: 

R2: 

2c 

S2 : 

U 

RC: 

R3 : 

S3: 

V 

RD: 

R4  : 

0M 

S4 : 

X 

RE: 

R5 : 

am 

S5 : 

Y 

RI: 

R6: 

'ms 

S6 : 

6xd 

R7 : 

°s 

S7 : 

A  A  ' 

m 

R8  : 

As 

S8: 

d 

R9 : 

ijSM 

S9 : 

f 

Initial  Flag  Status  and 

Use: 

0: 

OFF,  Unused 

2: 

OFF,  Unused 

1: 

OFF,  Unused 

3: 

OFF,  Unused 

Display  Status: 

DSP  4,  FIX,  DEG. 

User 

Control  Keys: 

A: 

f  AM 

a : 

Station  ID 

B: 

b : 

C: 

^S  *  A  s 

At 

D: 

d : 

m 


e : 


E:  Prepare  data  card 


Card  2 


Registers: 

RO: 

+  IDX 

SO: 

±  ID2 

RA: 

Rl: 

Atl 

SI: 

At2 

RB: 

R2: 

2cl 

S2 : 

2C2 

RC: 

R3 : 

S3: 

RD: 

R4 : 

6A1 

S4 : 

6A2 

RE: 

R5: 

XA1 

S5 : 

XA2 

RI: 

R6 : 

^Al 

S6 : 

^A2 

R7 : 

6B1 

S7 : 

6B2 

R8 : 

XBl 

S8 : 

XB2 

R9 : 

^Bl 

S9 : 

^B2 

Initial  Flag  Status  and  Use: 

0: 

OFF,  Vertex 

determination 

2:  OFF, 

Validity 

1: 

OFF,  Unused 

3:  OFF, 

Unused 

Display  Status; 

DSP  2,  FIX,  DEG 


User  Control  Keys; 


A:  a: 

B:  b: 

C:  c: 

D:  d: 

E:  Run  e: 


=  21295.87 
P 

= 1/298.26 


checking 
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Card  3. 


Registers : 

R0: 

ID1 

SO: 

ID2 

RA: 

M 

Rl: 

Atl 

SI: 

At2 

RB: 

u. 

R2 : 

(2c)  1 

S2: 

( 2c)  2 

RC: 

D, 

R3 : 

Al'  c  2 ' 

P 

S3: 

A2 

RD: 

Ao 

R4: 

0 ,  =0 

1  F 

S4 : 

62 

RE: 

a 

P 

R5 : 

>> 

H-* 

II 

*3 

S5: 

A2 

RF: 

f  = 

R6 : 

4 

S6  : 

^  2 

R7 : 

ci'  ci' 

H 

S7 : 

C2 

R8 : 

B.^  S/a 

=  r 

S8 : 

B2 

R9 : 

U1 

S9 : 

f  =  1/298.26 


Initial  Flag  Status  and  Use: 

0:  OFF,  Soln  A,  Soln  B  2:  OFF,  M/S  Vertex  Flag 

1:  OFF,  Unused  3:  OFF,  Unused 

Display  Status; 

DSP  2,  FIX,  DEG 


User  Control  Keys: 
A:  T, 


1 


B: 

C: 


a:  Soln  A 

b:  Soln  B 

c : 


639  *LBLE  21  13  Main  Routine:  Renewal  Solution 


1T>  '  O  '  O 

ro  ro  ro 


-r  co  in 


io  ^x>  'o 
ro  ro  ro 


j  in  o>  r  j  *u  >f  oj  to  ro  •*-*  to  10  1  u  ’h  r-i  io  •-«  ro  lo  o  •-« 

j  TJ.  i«J  t  in  Vi  V  t)  'T  111  vf  *t*  ’t  Irt  fo  *■•  V  lo  Q  Csl 

lit  I  *  I  * 


IO  lD  *0 

ro  — «  h) 


r>- 

-p  -J  p 
o  VJi 
ct  i>: 


o  r  -  *■'» 

•i*  o  — 4  — i  i 


Q  00  ID 
•I*  O  1  -P 
►-  O  CP 
co  cc  u: 


*-*  w  Q  to  »>j  o  Z  1 
o  u  j  o  :<  .j  n  : 
K*  Cl  U  U  U  O) 

CO  u;  Cic 


•-«  2:  w 

i  ^ 

o  to 

oc 


o  «- 

>  O 


•3-.  — .  c-j  r<»  .4*  m  to  r  ~  ..o  o-i  ca  — «  cm  ro  ^  in  ms  r  .  co  a,  co  <'i  ro  o-  i n  •-£.  r-  to  a,  c&  ro  Jp 
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D.  Loran-C  Fixing  Algorithms 

The  development  of  the  Loran  fixing  algorithms  in  this 
report  is  presented  in  more  detail  in  a  companion  report 
[Ref.  3]  and  will  not  be  repeated  here. 

The  basic  Loran-C  equation  [Ref.  4]  can  be  written  as 

T  =  [Tg  +  p(Tg)]  -  [Tm  +  p(TM)]  +  [Tb  +  p(TN)]  +  6  (1) 

where 

T  is  the  "indicated  time  difference"  in  microseconds, 
TM'TS  ;*'s  t*le  distance,  in  microseconds ,  from  the  master 
and  the  slave  to  the  receiver,  respectively, 

Tg  is  the  distance,  in  microseconds,  between  the 
master  and  the  slave, 

6  is  the  assigned  coding  delay,  in  microseconds,  and 
p (T)  is  the  secondary  phase  correction,  in  microseconds, 
for  an  all  sea  water  path  of  length  T. 

The  quantity 

At  =  [TB  +  P(TB) ]  +  6 

is  a  constant  for  each  master/slave  pair.  The  following 
World  Geodetic  System  1972  (WGS  72)  values  have  been  adopted 
for  Loran-C  navigation  [Ref.  4]: 
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Vq  =  299792458  me ter s/second  is  the  velocity  of  light 
in  free  space, 

n  =  1.000338  is  the  index  of  refraction  of  the  surface 
of  the  earth  for  standard  atmosphere  and  100kHz 
electromagentic  waves, 

ag  =  6378135.00  meters  is  the  equatorial  radius  of  the 
earth,  and 

f  =  1/298.26  is  the  flattening  factor  (1  -  b/ae,  where 
b  is  the  polar  radius)  of  the  earth. 

Accurate  formulas  for  computing  the  secondary  phase 
correction  p(T)  are  contained  in  Reference  4,  but  for  use 
with  the  handheld  calculator  the  following  linear  approximation 
[Ref.  3]  will  be  used: 


p  (T)  =  a.  +  a_T 


2  ' 


where 


and 


ax  =  -0.321, 
a2  =  0.000635. 


Using  this  approximation,  it  is  possible  to  solve  Equation  1 
for  the  quantity  Tg  -  T^.  We  find 


Ts  "  tm  =  (T  -  At) / ( 1  +  a2)  . 


(2) 


On  the  surface  of  a  sphere  a  hyperbolic  line  of  position 
can  be  represented  by  the  equation  [Ref.  3,  page  175] 
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J 


tan  r 


(3) 


_ cos  2a  -  cos  2c 

sin  2c  cos  w  +  c,  sin  2a 

where  the  origin  of  the  coordinate  system  is  at  the  prime  focus 
of  the  spherical  hyperbola,  2c  is  the  spherical  arc  joining 
the  foci,  2a  is  a  constant  for  any  one  LOP,  and  r  and  to 
are  the  spherical  coordinates  of  a  point  on  the  LOP.  If  the 
base  line  of  the  coordinate  system  is  the  arc  joining  the  foci 
then  to  is  the  spherical  polar  angle  from  the  base  line  to  a 
point  P  on  the  LOP  and  r  is  the  spherical  polar  distance 
(or  arc)  from  the  prime  focus  to  P.  Using  the  Loran  system 
we  take  c,  =  +1  if  the  prime  focus  is  at  a  master  station 
and  we  take  i,  =  -1  if  the  prime  focus  is  at  a  slave  station. 

If  we  take  v  =  vQ/n  to  be  the  velocity  of  100kHz  electro¬ 
magnetic  radiation  of  the  earth's  surface  then 


or,  using  Eq.  (2) , 


where 


2a  =  v(Ts  -  TM)/ae  , 


2a  =  (T  -  At) /a  , 


a  (i  +  a_) 

a  =  - — £_  =  21295.87  us  . 

P  v./n 


(4) 


The  baseline  between  master  and  slave  can  be  obtained  from 


2c  =  v  T  /a  .  (5) 

B  e 

Here  2c  is  computed  by  program  card  1  (preparation  of  master 
data  cards)  using  the  algorithm  in  Section  E. 
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Consider  a  Loran-C  triplet  with  master  stations  colocated. 


Let  and  £2  denote  the  azimuth  angles  of  slave  1  (S^ 

and  slave  2  (S2) ,  respectively,  measured  from  North  toward  the 
East  from  the  master  stations  (M)  (see  Fig.  3) .  Further,  let 
a  and  r  denote  the  azimuth  and  spherical  polar  arc  (distance) 
of  the  receiver  (R)  from  M.  For  this  geometry,  Eq.  (3)  can 
be  written  as 


tan  ri  ”  C~~cosTo^T~5  +~a7 


where 


A.  =  <; .  sin  2a. 
i  i  i 


B.  =  cos  2a.  -  cos  2c. 
ill 


C .  =  sin  2c . 
l  l 


for  the  i—  Loran  pair,  i  =  1,2.  Since  r  =  ^  =  r2,  tan  r^ 
can  be  eliminated  in  Eq.  (6) .  The  resulting  equation  can  be 


rewritten  as 


where 


C  cos  a  +  S  sin  a  =  K  , 


C  =  BlC2  cos  C2  -  B^  cos 
S  =  B^C2  sin  S2  -  B.C,  sin 
k  =  b2a^  -  b^a2  . 


If  we  define  p  >  0  and  y  by  the  equations 


then 


P  cos  y  =  C, 
p  sin  y  =  S, 


P  =  /c2  +  S2  , 


y  =  qatn(S,C) 


Here  the  function  qatn(y,x)  is  the  arctangent  of  y/x 
adjusted  for  the  proper  quadrant  according  to  the  signs  of  x 
and  y.  A  compact  form  of  this  function  is 


qatn(y,x)  =  tan 


x  +  10  t(x  =  0?) 


+  TTt(x  <  0?) 


where 


t(z)  =  1  when  z  is  true 


t ( z)  =0  when  z  is  false. 


When  convenient  we  will  use  the  notation  qatn(y/x)  inter¬ 
changeably  with  qatn(y,x).  The  qatn  function  is  equivalent 
to  the  polar  angle  obtained  using  the  rectangular  to  polar 
conversion  function  on  the  HP-67/97. 

Now  substitute  Eq .  (8)  into  Eq.  (7)  and  solve  for 


ot  =  y  +  cos  (k/p  ) 


to  obtain  the  azimuth  angle  a  of  the  two  points  of  intersection 
of  the  LOP's.  Finally  we  obtain  a  value  for  r  by  substituting 
each  a  into  Eq.  (5) .  We  find  that 


r  =  qatn 


cos  (a-E;^)  + 


for  i  =  1  or  2 . 


The  distance  and  azimuth  from  M  or  the  triplet  vertex  can  be 
converted  into  the  latitude  and  longitude  of  the  two  possible 
positions  of  R. 


The  fixing  algorithm  then  uses  a  and  r  in  the  d 
solution  algorithm  of  spheroidal  geodesy  (Section  F) 


trert 


The  Reverse  (Inverse)  Solution  Algorithm 


This  reverse  solution  algorithm  is  a  modification  of  the 
first  order  in  flattening  (f)  algorithm  given  by  Thomas  [Ref.  5, 
pp.  8-10].  Thomas'  notation  has  been  followed  as  closely  as 
possible  for  ease  of  comparison  of  the  algorithms.  The  qatn 
function  is  defined  in  Section  D.  West  longitudes  (X)  and 
South  latitudes  (!)  are  negative.  We  are  given  the  points 
Pl^  l'^l^'  on  t*ie  sPher°id  and  are  to  find  the 

distance  S  between  the  points  and  the  forward  and  back 
azimuths,  and  l21*  Given  quantities  are  cf  ^ ,  X^,  i}i 

and  *  .  No  assumptions  about  the  relative  location  of 

and  P  are  required.  The  modified  -vmv  solution  algorithm 


-  tan  [(1-f)  tan  $  .  )  ,  i  =  1,  2, 

,  +  •<  J  /2,  An  =  («.-•• ,  )/2.  ,M  =  \  -  \  , 

m  1  2  m  2  I  2  1 

2  j  2  2 

r  2  ■  '2,  H  =  i'i) s  -  sm"  h  =  cos"  •*  -  sin"  At)  =  cos  0  cos  0  , 

m  m  m  m  m  1  2 


I.  -  sin  '0  +  H  sin  A' 

m  rn 


;in"(d/2),  1  -  L  =  cos"  (d/2), 


i  -  cos  (1  -  21,),  U  =  2  sin  II  cos  2  0  /(I  -  L)  , 

tn  m 


sm  A‘<  cos  -  -I,,  X  =  ll  ♦  V,  V  -  U  -  V, 
m  m 


I  i  in  d,  ■  d  -  t  (TX  -  Y)/4,  S  =  a  ( T-c  . cJ )  sin  d, 


i'  2  (Y  -  (  1  -  21.)  (4  -  X)  I  ,  A  =  fT/2, 


t,  =  qatn(-sin  AS  cos  AX',  cos  6  sin  AA')f 
1  m  m  m  tn 


t„  =  qatn(cos  A0  cos  AX',  sin  0  sin  AX')  , 
2  m  m  m  m 


a!2  =  fcl  +  V  a21  =  fcl  '  fc2 


This  reverse  solution  algorithm  is  used  by  program  card  1 
(preparation  of  master  data  cards)  to  compute  the  baseline 
distance  2c  and  the  azimuths  £MS  and  between  the 

master  and  slave  stations  of  a  Loran  pair. 

Details  of  the  modifications  made  to  Thomas'  algorithm 
are  contained  in  Reference  3. 
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F.  The  Direct  Solution  Algorithm 


This  direct  solution  algorithm  is  a  modification  of 
the  first  order  in  flattening  (f)  algorithm  given  by  Thomas 
[Ref.  5,  pp.  7-8].  Thomas'  notation  has  been  followed  as 
closely  as  possible  for  ease  of  comparison  of  the  algorithms. 
The  qatn  function  is  defined  in  Section  D.  West  longitudes 

and  South  latitudes  are  negative.  We  are  given  the  point 

Pl^l'*l*  °n  the  sP*herio;i-d'  where  *  are  the  geodetic 

latitude  and  longtiude  (geographic  coordinates) ;  the  forward 
azimuth  and  the  distance  S  to  a  second  point  P2^2'*2^ 

and  from  these  we  are  to  find  the  geographic  coordinates 
^2#  Ap  and  fche  back  azimuth  a21'  T^e  9aven  quantities  are 
■t  ,  A  ,  1^2  and  S.  assumptions  about  the  relative  loca¬ 
tion  of  and  P2  are  required.  The  modified  direct 

solution  algorithm  is: 

=  tan  P  [  <  1  —  f )  tan  v  ^ !  ,  M  =  cos  6^  sin  « 

N  =  cos  < * ^  cos  u  c ^  =  fM,  c ^  =  f ( 1  -  M2)/4, 

D  =  1  -  2c  ^  ~  c^M,  P  =  C'^/D ,  0  ^  =  <;atn(N,  sin  6  ) 


d  = 

S/ ( a  D) , 

e 

u  : 

=  2(a 

Yd) 

,  V7  = 

1  -  2P 

cos  u , 

V  = 

cos ( u  +  d ) , 

Y  = 

2PVW 

sin  d , 

Act  = 

d-Y , 

i  - 

21 

;ntn (-M,  - 

(N 

cos 

A  _ 

sin  0 

sin  An )  ]  , 

K  = 

(1-f) (M2  + 

(N 

cos 

A  •  - 

sin  (‘ 

sin  A  a 

)2i1/2  , 

> 

tan  l  ( s l n 

t  i 

1 

cos 

An  + 

N  sin 

)/K] 

t 

qnt  n  (  s  i  n  ' 

s 

in  •  i 

\  2  '  r  <’S  ]  •  os  Ai"  * 

sin  (i  sin 

H 

1  '  • 

'•/I  - 

II, 

1  +  "V* 

\ 


This  direct  solution  algorithm  is  used  by  program 
card  3  (improved  fix  program)  to  compute  the  latitude  and 
longitude  of  the  receiver  using  the  azimuth  and  range  of  the 
receiver  from  the  Loran  triplet  vertex. 

Details  of  the  modif ications  made  to  Thomas'  algorithm 
are  contained  in  Reference  3. 


G. 


Discussion  and  Some  Typical  Results 


f 


The  HP-67  program  design  specifications  of 
COMPATWINGSPAC  [Ref.  6]  are  contained  in  the  following 
statement . 


"There  is  a  need  for  an  HP-67  program  that  will  compute  a 
geographical  position  from  two  Loran  delay  rate  readings. 

Several  methodologies  are  available  to  compute  the 
desired  position  but  computational  complexities  increase 
with  the  desired  accuracy  and  flexibility.  The  most 
desirable  accuracy  would  be  an  error  of  less  than  4  n.mi. 
at  a  range  of  500  n.mi.  wit,,  less  error  closer  to  the 
stations.  It  is  likely  that  program  length  considerations 
will  require  that  the  station  pairs  have  a  common  site 
(i.e.  two  slaves  or  two  masters  at  the  same  location). 

This  is  not  an  unusual  situation  as  evidenced  by  strinqs 
of  station  pairs  along  coast  lines.  A  data  card  will 
probably  be  necessary  for  the  station  pairs  to  be  used. 

However,  more  than  one  proqram  card  is  unacceptable  due 
to  the  decrease  in  functional  utility  when  compared  to 
the  manual  plottinq  method.  As  a  final  requirement,  the 
fix  should  be  obtainable  on  either  side  of  the  baselines 
connecting  the  stations,  and  not  limited  to  a  geometric 
position  relative  to  one  side  or  the  other  of  the 
stations . " 

It  was  further  stated  that  the  maximum  computation  time  to 
obtain  a  fix  be  1.5  minutes. 

It  is  felt  that  these  design  goals  have  been  satisfied. 
Although  one  program  is  required  to  prepare  master  data  cards 
for  all  Loran-C  pairs  and  a  second  card  is  required  to  prepare 
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operational  data  cards,  one  each  for  every  triplet,  this  pre¬ 
paration  should  be  done  only  once.  The  data  cards  should  be 
supplied  to  users  verified  and  labeled,  by  the  Fleet  Mission 
Program  Library.  One  program  card  and  an  appropriate  oper¬ 
ational  data  card  are  all  that  is  required  for  the  fixing 
algorithm. 

The  fixing  algorithm  will  display  one  of  the  two 
possible  receiver  positions  in  38  seconds  following  the  entry 
of  the  time  delay  readings.  Since  there  are  situations  in  which 
either •  of  the  two  solutions  could  be  the  valid  solution;  the 
decision  of  which  solution  to  use  should  be  left  to  the  operator, 
not  the  program  designer. 

Testing  of  the  algorithm  for  all  Loran-C  triplets  and 
positions  relative  to  those  triplets  was  too  extensive  a 
program  to  be  carried  out  in  the  available  time.  Some  "typical" 
scenarios  however  are  presented  in  Tables  I  through  IV.  As  can 
be  seen  all  errors  are  all  well  within  the  design  specif ications 
of  4  n.mi  at  500  n.mi  range  from  the  stations.  The  time  delay 
values  in  these  Tables  were  generated  using  a  program  discussed 
in  Reference  3.  It  is  recommended  that  the  P-3  community  test 
the  algorithm  for  accuracy  in  known  areas  of  operation  and 
examine  the  results  for  possible  regions  in  which  the  algorithm 
may  fall  outside  the  design  requirements.  Such  testing  should 
be  compatible  with  the  known  "unreliable  regions"  shown  on  the 


Loran-C  charts. 


Table  I.  Moffett  Field  South 


Position 

Lat  Long 

Indicated 

9940X 

Time  Delay 

9940Y 

Lat (N) 

Fix 

Long  (W) 

Error 

n.mi 

24°N 

122  °W 

27726. 19 

40912.76 

23°59 

55" 

122°00‘01" 

0.08 

26 

122 

27715.97 

40998. 39 

25°59 

57" 

122  °00 ’01" 

0.05 

28 

122 

27702.41 

41117.84 

27°59 

59" 

122°00 ' 00" 

0.02 

30 

122 

27683.53 

41291.85 

29°59 

59" 

122  °00 1 00" 

0.02 

32 

122 

27655.47 

41555.46 

32  °00 

00" 

122"00 ' 00" 

0.00 

34 

122 

27609.63 

41959.57 

34  °00 

00" 

122  °00 ' 00" 

0.00 

36 

122 

27523.56 

42544.11 

36°00 

00" 

121°59' 59" 

0.01 

38 

122 

27334.61 

43248.22 

38°00 

00" 

121°59' 58" 

0.03 

! 


Table  II.  Moffett  Field  West 


- — - - — 

Position 

Lat  Long 

Indicated 

9940Y 

Time  Delay 

9940W 

Lat (N) 

Fix 

Long (W) 

Error 

n.mi 

37°N 

122  °  W 

42892.86 

16257.23 

36°59 

59" 

122 “00 '01" 

0.02 

37 

125 

43056.68 

15765.13 

37°00 

00" 

125  °00 1 00" 

0.00 

37 

128 

43137.78 

15327.12 

37°00 

00" 

128°00 ' 00" 

0.00 

37 

131 

43191.10 

14970. 77 

37°00 

00" 

1 31°00 ' 00" 

0.00 

37 

134 

43232.38 

14683.74 

37  °00 

00" 

1 34  °00 ' 00" 

0.00 

37 

137 

43267.42 

14449.40 

37°00 

00" 

1 37°00 ' 00" 

0.00 

37 

140 

43298.80 

14254.02 

37  °00 

00" 

140°00 ' 01" 

0.01 

37 

143 

43327.85 

14087.43 

37  °00 

01" 

142°59' 59" 

0.02 

Table  III.  Brunswick  Northeast 


Position 

Lat  Long 

Indicated 

7930Z 

Time  Delay 

9930X 

Lat (N) 

Fix 

Long (W) 

1  Error 

n.mi 

60  °N 

30°W 

52437.86 

28451. 72 

60°00 ' 03" 

29°59' 32" 

0.24 

58 

35 

51960.93 

28391. 50 

o 

o 

o 

o 

o 

CO 

in 

34059' 46" 

0.11 

56 

40 

50992.37 

28359. 15 

55°59 1 59" 

39°59 ' 54" 

0.06 

54 

45 

49292.46 

28370. 85 

53°59 ‘ 59" 

44°59' 57" 

0.03 

52 

50 

47165.60 

28490. 64 

52°00 1 00" 

49°59 ' 59" 

0.01 

50 

55 

45236.59 

29070. 48 

50°00 ' 00" 

55°00'00" 

0.00 

48 

60 

44505.60 

30991. 94 

48°00  *  00" 

60°00 ' 00" 

0.00 

46 

65 

44475.70 

l 

33697. 14 

46600'00" 

65°00'00" 

0.00 

44 

70 

44588.91 

36567. 42 

43°59' 59" 

69  °59' 59" 

0.02 

Position 


Table  IV.  Jacksonville  Southeast 

Indicated  Time  Delay  F 


Error 


Lat 

Long 

9930W 

9930X 

Lat(N) 

j  Long(W) 

n.mi 

9°N 

47°W 

13058.04 

36466.46 

8°59 ' 19" 

46°59 

’22" 

0.92 

12 

52 

12984.71 

37288.  35 

11°59' 34" 

51°59 

’37" 

0.57 

15 

57 

12898.73 

38267. 58 

14°59 ' 44" 

56°59 1 

'47" 

0.34 

18 

62 

12793.91 

39431.  32 

17°  59 ' 52" 

61°59 1 

’54" 

0.16 

21 

67 

12656.52 

40794.  36 

20°59 ' 56" 

66°59 1 

1  57" 

0.08 

24 

72 

12451.30 

42330. 55 

23°59' 59" 

71°  59 1 

’59" 

0.02 

27 

77 

12097.12 

43876.  62 

27°00 ' 01 " 

77°00 1 

00" 

0.02 

30 

82 

12973.95 

44768. 53 

30°00’01" 

82  °00 1 

06" 

0.09 
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APPENDIX.  Loran-C  Station  Parameters 

The  following  list  contains  the  pertinent  parameters 
for  each  Loran-C  station  pair.  This  list  was  compiled  from 
data  in  Reference  4 .  Each  column  contains  the  following 
information: 

1.  The  Loran-C  station  pair  designator 

2.  At,  the  sum  of  the  coding  delay  plus  one  way  baseline 
time,  including  the  secondary  phase  correction  for  an 
all  seawater  path,  in  microseconds. 

3.  The  master  station  latitude. 

4.  The  master  station  longitude. 

5.  The  slave  station  latitude. 

6.  The  slave  station  longitude. 

In  this  list,  negative  longitudes  are  West  and  positive 

longitudes  are  East.  If  desired,  this  convention  may  be 

reversed  since  the  algorithms  are  independent  of  such  external 

conventions;  if  this  is  done,  care  should  be  taken  that  the 

signs  of  all  longitudes  in  the  list  are  reversed.  In  columns 

3  through  6  the  latitudes  and  longitudes  appear  to  be  in 

decimal  form,  but  the  actual  format  is  DDD.MMSSFF  (which  is 

compatible  with  the  HP-67/97  H.MS  input  mode)  where 

DDD  designates  degrees, 

MM  designates  minutes, 

SS  designates  seconds,  and 

FF  designates  hundredths  of  seconds. 
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4990X 

15972. 

ul  o'  .* 

16. 

444395, 

-163. 

303128, 

20. 

144916 

4990V 

34253. 

18, 

16. 

444395, 

-169. 

303120, 

28. 

234177 

5930X 

13131. 

88 

46. 

482720, 

-067. 

553771, 

41. 

151133 

5930V 

28755. 

02 

46. 

482720, 

-067. 

553771 , 

46. 

463218. 

5990X 

13343. 

60 

51. 

575878, 

-122. 

220224, 

55. 

262085 

5990V 

28927. 

36 

51. 

575878, 

-122. 

220224, 

47. 

034799 

5990Z 

42266. 

6-3 

51. 

5  c 58 r 8 , 

-122. 

220224, 

50. 

362972 

7930W 

15068. 

02 

59. 

591727, 

-045. 

102747, 

64 . 

542658 

7930X 

27803. 

<  •** 

59. 

591727, 

-045. 

102747, 

62. 

175368 

7930Z 

48212. 

20 

59. 

591727, 

-045. 

102747 

46. 

463218 

7960X 

1 3804 . 

45 

63 . 

194281, 

-142. 

4831 90 

C7 
■J  l  . 

262021 

7968V 

2965 1 . 

14 

63. 

194281, 

-142. 

483190 

55. 

262085 

7970W 

30965. 

64 

62. 

175968, 

-007. 

04267 1 

54. 

482980 

7970X 

15048. 

10 

62. 

175968, 

-007  . 

042671 

68 . 

3806 1 5 

7978V 

48944. 

cr  - 

62. 

175968, 

-007. 

04267 1 

64 . 

542658 

7970Z 

63216. 

30 

62 . 

175963, 

-007. 

84267 1 

70. 

545261 

7980W 

12309. 

54 

30 . 

593874 , 

-085. 

1 00330 

30 . 

433302 

7988X 

27*42. 

38 

30. 

53-3874  • 

-085. 

1 009:90 

26 . 

315501 

7980V 

45201 . 

ij 

30 . 

593874 , 

-085. 

1 00930 

2 1  • 

015849 

7938Z 

61542. 

"? 

30. 

593874 , 

-085 . 

100330 

34. 

034604 

7 999X i 

i  -7 err 

J.  2.  i  -J . 

97 

j"' 

522061 . 

016. 

430596 

35. 

312088 

7  9  90V 

Z-  2  Cm  1  Z'  * 

30 

i~i 

522961 , 

016. 

430596 

46. 

532095 

79907 . 

50?  9*3 . 

69 

o 

52206 1 , 

016. 

430596 

42’  • 

033649 

•:;978K 

14355. 

i  i 

510754, 

-087 . 

291214 

30. 

593874 

970!  •; 

31162. 

06 

39. 

510754, 

-087. 

291214 

42. 

42506© 

3970’’ ' 

47753. 

74 

3 '9 . 

510754, 

-687 . 

291214 

48. 

364934 

*C'0t: 

t‘605. 

51 

34. 

034604, 

-077. 

544676 

-’7 

015849 

7i  7  '  ' 

'  i  "*  m 

66 

34. 

934604 ■ 

-077. 

544676 

46 ! 

4 632 1 8 

7; 

52541. 

•~i  i 

34. 

034604 , 

-077. 

544676 

41. 

151193 

”  -7*  “•  t~, 

63560. 

7  •*' 

34. 

034604 , 

-077. 

544676 

99 

516754 

'  '49;j. 

1 3796. 

_-i  0 

39 . 

330662 , 

-110. 

495637 

47. 

034799 

’  ?4i'' . 

7  '*!0'4  4 

50 

99 . 

9.30662  ■ 

-113. 

495637 

465699 

“i  “  ^  •  7(1  i 

■  4  i  ?-:"7 , 

30 

39 . 

330662 , 

-118. 

495637 

25. 

191818 

7.  *.  » 7(1 

!  I-."7 ?7. 

28 

42.’ 

425060 , 

-076. 

493386 

46. 

482720 

■'960;, 

3  f7  *“"•  3 . 

9  3 

42. 

425060, 

-076 . 

493386 

41. 

151193 

9  :'6  j".' , 

-  * "  ‘ *'  1 . 

65 

42. 

425060 ■ 

-076 . 

493386 

34. 

034604 

99t  8Z 

cr  —  i  r 
»  J.  -  -  c~  • 

0  C 

42. 

425060, 

-076. 

493386 

39. 

510754 

15233. 

94 

24. 

4304 1  , 

141. 

19290 

24. 

17077 

7  f  rz  ■“'«  cr 

12 

24. 

48041  .• 

141. 

1 9290 

42. 

443700 

”  '4  7  pi 

■  5? 433. 

4  Z< 

1  ^ 

24. 

4804 1  , 

141. 

1 9290 

26. 

362499 

“970  Z.. 

397<£. 

“■*  t| 

24. 

48041  , 

141. 

i  3290 

09. 

324566 

:,Qp('  ■ 

14075. 

. ’i . 

090988 • 

- 1 70 . 

145981 

52. 

494505 

• -V 

;'.,069. 

0C 

cr-» 

»  ■ 

090988 , 

-170. 

145981 

65. 

144912 

.-j  *  “if---* 

46598. 

10 

cr 

■_» »  , 

09W988 , 

-170. 

145981 

57. 

262021 

-155.530976 
-178. 

-069. 

-053. 

-131. 

-119. 

-127. 

-023. 

-007. 

-053. 

-152. 

-131. 

+008. 

+014. 

-023. 

-008. 

-090. 

-097. 

-080. 

-077. 

012. 

027. 

003. 

-085. 

-076. 

-094. 

-080. 

-053. 

-069. 

-087. 

-119. 

-122. 

-114. 

-067. 

-069. 

—0  ( ( . 

-087. 

153. 

143. 

128. 

138. 

+  173. 

-166. 

-152. 
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Station 


Location 


4990 

Central  Pacific 

5930 

East  Coast,  Canada 

5990 

West  Coast,  Canada 

7930 

North  Atlantic 

7960 

Gulf  of  Alaska 

7970 

Norwegian  Sea 

7980 

Southeast  U.S.A. 

7990 

Mediterranean  Sea 

8970 

Great  Lakes 

9930 

East  Coast,  U.S.A. 

9940 

West  Coast,  U.S.A. 

9960 

Northeast  U.S.A. 

9970 

Northwest  Pacific 

9990 


North  Pacific 
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New  London,  CT  06320 

Naval  Ship  Engineering  Center 
Hyatts vi lie,  MD  20782 

Naval  Coastal  Systems  Laboratory 
Panama  City,  FL  32401 

Naval  Air  Systems  Command 
Code  370 

Washington,  D.C.  20361 

Naval  Sea  Systems  Command 
Code  0  3424 

Washington,  D.C.  20362 

Naval  Underwater  Systems  Center 
Newport,  Rhode  Is  la  id  02  84  0 

Naval  Ordnance  Station 
Indian  Head,  MD  20140 

Naval  Surface  Weapo  is  Center 
Dahlgren,  VA  22448 

Anti-Submarine  Warfare  Systems  Project  Office 
Department  of  the  Navy 
Washington,  D.C.  20360 
Attn:  Code  ASW-137 

Surface  Warfare  Development  Group 
Naval  Amoh ibious  Base 
Little  Creek,  VA  23511 
Attn:  Code  N32  and  C.  Cartledge 

COMMA  VS  U  RF I  .ANT 
Norfolk,  FA  23511 
Attn:  Code  N’t  12 

Deputy  Cerumen der 

Opcr at  n ml  Ti  nt  and  Evaluation  Force  Pacific 

NAS  North  Island 
San  Diego,  OA  92135 
Attn:  Code  202 

VT-24 

NAS 

Jacksonville,  FL  32212 
Attn:  i. TJO  Kenneth  D.  Walker 
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Assistant  Wing  TAC  D&E  Officer  1 

ASW  Operations  Department 

NAS 

Cecil  Field,  FL  32215 
Attn:  Code  VS-24 

Naval  Electromagnetic  Spectrum  Center  1 

Naval  Communications  Unit 
Washington,  D.C.  20390 
Attn:  CDR  Claude  I aVurre 

VS-Support  Unit  1 

NAS 

Cecil  Field,  FL  32215 
Attn:  CWO  Wentworth 

P ROC AL- Code  18  1 

Fleet  Readiness  Office 

NOSC 

San  Dicqo,  CA  92152 
Attn:  Jim  Grant 

Chief  of  Naval  Materiel  1 

MAT  0  81)1 

Department  of  the  Navy 
Washington,  D.C.  20360 
Attn:  CDR  Philip  Harvey 

LT  Jan  Smith  3 

COMCRUDESGRU  TWO 

FPO  New  York,  NY  09501 

Dr.  Thomas  Burnett  1 

Dept,  of  Manaqemont  Science 

OR  Study  Group 

U.S.  Naval  Academy 

Annapolis,  MD  21401 

U.  S.  Naval  Oceanoqraph ic  Office  1 

Washington,  D.C.  20309 
Attn:  Paul  D.  Thomas 

Navigation  Department  1 

Defense  Mapping  Aqency 

Hydrograph ic  'Topograph ic  Center 

Washington,  D.C.  20315 

Attn:  G.  R.  Deyoung 
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LT  Michael  D.  Redshaw 

Helicopter  Antisubmarine  Squadron  7 

FPO  New  York  09501 

Mr.  Raymond  F.  Fish 

Naval  Underwater  Systems  Center 

Newport,  RI  02840 

LT  Kenneth  W.  Peters 
Commanding  Officer 
U.  S.  Naval  Air  Facility 
APO  New  York  09406 

LT  Peter  W.  Marzluff 

Surface  Warfare  Officer's  School 

Newport,  RI  02840 

LT  Ray  C.  Pilcher,  Jr. 

Surface  Warfare  Officer's  School 
Newport,  RI  02840 

LT  William  M.  Yerkes 
ASWOC  NAF 
S’gonella  Sicily 
FPO  New  York  09523 

LT  A.  J.  Koorey 
Commander 

Cruiser-Destroyer  Group  Two 
FPO  New  York  09501 

LCDR  E.  G.  Schwier 
USS  Estocin  ( FFG  15) 

FPO  New  York  09501 

LC DP  R.  R.  Hi  1  Iyer 
Code  N311 

Commander,  Carrier  Group  Four 
FPO  New  York  09501 

LCDR  J.  W.  Pattison 

Code  N32  1 

Commander,  Carrier  Group  Eight 
FPO  New  York  09501 

LT  Richard  C.  Myers 
ASWOC  NAF  LAJES 

A cores 

APO  New  York  09523 
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LT  R.  N.  Christianson 
ASW 

Naval  Postgraduate  School 
Monterey,  CA  93940 

LT  L.  R.  Erazo 
ASWOC  SIGONELLA 
FPO  New  York  09523 

LT  R.  E.  Springman 
USS  Lexington  (ANT-16) 

Naval  Air  Station 
Pensacola,  FL  32504 

LT  Philip  D.  Ward 
Patrol  Squadron  26 
Special  Project  Det 
Naval  Air  Station 
Brunswick,  ME  04011 

LT  Robert J.  Knight 

USS  John  F.  Kennedy  (CV-67) 

FPO  New  York  09501 

LT  Richard  0.  White,  Jr. 

USS  Dwight  D.  Eisenhower  (CVN-69) 
FPO  New  York  09501 

LT  David  E.  Hebdon 

USS  John  F.  Kennedy  (CU-67) 

FPO  New  York  09501 

LCDR  Michael  A.  Santoro 
Air  Antisubron  41 
NAS  North  Island 
San  Diego,  CA  92135 

LCDR  Robert  M.  Hanson 
Helasron  Four 
NAS  North  Island 
San  Diego,  CA  92135 

LT  Michael  D.  Thomas 
2476  Ridgecrest  Ave . 

Orange  Park,  FL  32073 

LT  Carl  E.  Garrett,  Jr. 

Class  65,  Dept.  Head  Course 
Surface  Warfare  School,  Naval  Base 
Newport,  RI  02840 
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Naval  Postgraduate  School 
Monterey,  CA  93940 

Attn:  Prof.  R.  N.  Forrest,  Code  55Fo 

Prof.  R.  H.  Shudde,  Code  55Su 
R.  J.  Stampfel,  Code  55 


LT  M.  D.  Clary 
USS  PAIUTE  ( ATF-159 ) 
FPO  New  York  09501 
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